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Abstract 

Family symmetry could explain large mixing of the atmospheric neutrinos. The same symmetry 
could explain why the flavor changing current processes in supersymmetric standard models can 
be so suppressed. It also may be able to explain why the proton is so stable. We investigate these 
questions in a supersymmetric, renormalizable extension of the standard model, which possess a 
family symmetry based on a binary dihedral group Qg • We find that the amplitude for /x — > e + 7 
enjoys a suppression factor proportional to|(yMNs)e3| — nie/ (\/2m^) ~ 3.4 x 10~^, and that B{p 
fi^)/ B{p K^e~^) ~ |(^MNs)e3p — 10~^, where Vmns is the neutrino mixing matrix. 
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I. INTRODUCTION 



The remarkable success of the standard model (SM) suggests that we have a highly non- 
trivial part of a more fundamental theory for elementary particle physics. In spite of this 
success, the SM suffers from various problems. One of them is that due to the quadratic 
divergence of the Higgs mass its natural scale can not exceed 0(1) TeV ^|. Therefore, 
in order to extend the SM in a natural way, the quadratic divergence has to be cancelled 
As it is well known today, low energy supersymmetry (SUSY) is introduced to protect 
the Higgs mass from the quadratic divergence 0,1^. Unfortunately, SUSY is broken, and 
therefore its breaking should be soft to maintain the very nature of low energy SUSY, 
whatever its origin is 0, 0]- 

Renormalizability allows an introduction of a certain set of soft supersymmetry breaking 
(SSB) parameters. In the minimal supersymmetric standard model (MSSM), more than 
100 SSB parameters can be introduced j^. Moreover, these parameters should be highly 
fine tuned so that they do not induce unacceptably large flavor changing neutral currents 
(FCNCs) and CP violations 0, Q B S E3- Why should they be so fine tuned? What 
controls them ? This is the so-called SUSY flavor problem, which has existed ever since 
SUSY found phenomenological applications [ll|. 



There are several theoretical approaches to overcome this problem 0, 0|01"|2B- 
of them is to base on a family symmetry js^j-j^]. However, if the family symmetry is 



hardly broken (i.e. broken by operators with dimension > 4), non-symmetric SSB terms 
can be generated in higher orders in perturbation theory. Then it is not possible to make 
quantitative statements on the induced non-symmetric SSB terms, because these terms are 
infinite in perturbation theory, however small the hard breaking terms are. In this paper, 
we consider a supersymmetric model, recently proposed in [26|, which posses a nonabelian 
discrete family symmetry based on a binary dihedral group Qq. We investigate how the 
SUSY flavor problem can be softened by Qq family symmetry, which is our first task. We 
find that except for few cases the fine tuning of the SSB parameters is not necessary in this 
model. The only observable process would be /i ^ e -|- 7, whose branching fraction enjoys a 
suppression factor proportional to|(VMNs)e3p — fnl/i'^fnfi) — 10~^, and is three (four) orders 
of magnitude larger than that of r — e(/i) + 7 in the model, where Vmns is the neutrino 
mixing matrix. 

In the MSSM without R parity invariance, there exist about 100 dimension-four operators 
that violate the baryon number or lepton number conservations. Fortunately, R invariance 
can kill all these dangerous dimension-four operators, while it together with B — L conser- 
vation allows two types of dimension- five operators that lead to proton decay ji^, E3- 
is therefore widely expected that these operators are responsible for an observable proton 
decay. Indeed, in the minimal SUSY grand unified theory (GUT) based on SU{5) jlH . 
these operators force protons to decay faster than the experimental bound [2^. Given a 
family symmetry, the flavor structure is fixed, at least partially, and it can control the flavor 
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changing effects in proton decay, too '30]. In this paper we assume that some unknown 
Planck scale physics, which respects Qq family symmetry, induces dimension-five operators 
that lead to proton decay. We first would like to find out all B and L violating opera- 
tors with dimension < 5 allowed by Qq, and then to investigate how the family symmetry 
acts on proton decay. If we assume the degeneracy of the scalar quark masses, then the 
branching fraction B{p ^ e"*") is similar to that of ji^l that it becomes comparable 
with B{p ^ z/). If the degeneracy assumption is dropped, B(j> ^ e"*") can become 
0(10'^) X B{j) — 5> u), which is still two orders of magnitude lager than that of the 
minimal SUSY GUT case. It turns out that in contrast to the case of [301] the branching 
fraction B{p jX^), which turns out to be proportional to |(VMNs)e3p, is five orders of 

magnitude smaller than B{p e^) in the present case, irrespective of the degeneracy 

of the scalar quark masses. The main reason for this is the maximal mixing in the charged 
lepton sector, which serves the maximal mixing of the atmospheric neutrinos in this model. 
We also give an upper bound on the coupling constants that are assumed to be generated 
by Planck scale physics. We find that the family symmetry can raise the upper bound by 
two orders of magnitude. 

In sect, II we outlme the model of Q. Because of the uonabelian family symmetry, the 
Yukawa sector of the model is strongly constrained, and the redundancy of the parameters 
of this sector is significantly reduced. As a consequence, we can explicitly give the unitary 
matrices, Ul and Ur, that rotate the left and right-handed fermions to diagonalize the 
corresponding mass matrices. They are needed to quantitatively discuss the SUSY flavor 
problem in the super CKM basis in sect. III. For completeness, we re-present the predictions 
from Qq family symmetry. In sect. IV we consider the B and L violating operators that 
are allowed by the family symmetry and then calculate the dominant proton decay modes 
as function of superpartner masses. We conclude in sect. V. 



II. THE MODEL 

A. Group theory of Q2N 

The binary dihedral group Q2N {N = 2, 3 , . . . ) is a finite subgroup of SU{2) ^ . 
Its defining matrices are given by j26|, |3l| 

- l^-sin^ cos^j' 1^0 -zj' 

where On = 2tt/N . Then the set of AN elements of Q2N is given by 

Gq2N = {R2N-, {R2nT-, • • • , {R2nY^ = 1, R2nPq, {R2nTPq-, • • • , {R2Nf'^ Pq = Pq} ■ (2) 



^ Models based on dihedral flavor symmetries, ranging from Z33(~ 6*3) to Qq and Dj, have been recently 
discussed in 0, [ii-lii. 
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Only one- and two-dimensional irreps. exist. There exist four different one-dimensional 
irreps of Q2N1 which can be characterized according to Z2 x Z4 charge: 



l+,o, l-,o, l+,2, l-,2 for iV = 2,4,6,..., 

l+,o, l+,2, l-,3 for iV = 3,5,7,..., (3) 

where the l+,o is the true singlet of Q2N, and only 1_ 1 and 1_,3 are complex irreps. The 
— 1 different two-dimensional irreps are denoted by 

2,, £=!,.. .,iV-l. (4) 

2^ with odd £ is a pseudo real representation, while 2^ with even £ is a real representation. It 
is straightforward to calculate the Clebsch-Gordan coefficients for tensor products of irreps. 
The following multiplication rules in Qq in particular are used to construct the model |2q . 

l+,0 + l+,2 + 22 

-aihi + 0262 \ (5) 
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(6) 



B. Qq assignment and superpotential 

In Table H] we write the assignment of the quark, lepton and Higgs chiral supermul- 
tiplets ^, where Q, Qs, L, L3 and if", , if'', stand for SU{2)l doublets supermultiplets 
for quarks, leptons and Higgs bosons, respectively. Similarly, SU{2)l singlet supermultiplets 
for quarks, charged leptons and neutrinos are denoted by f/^, D'^, D3, and A^'^, N^. 

S,T and Y are SU{3)c x SU{2)l x U{1)y singlet Higgs supermultiplets. As an alterative 



to the abelian discrete R symmetry Z12R of [26[, for which one needs more singlets to con- 
struct a desired Higgs sector, we introduce an abelian Z4 symmetry. The Z4 constrains the 
Higgs sector, where it does not constrain anything in the Yukawa sector. When discussing 
proton decay in sect. IV , we will assume neither nor Zi2r., because the symmetry of the 
Higgs sector may have a stronger model dependence. We then write down the most general, 
renormalizable, Qq x Z4 x R invariant superpotential W: 

W = Wq + Wl + Wh, (7) 



The same model exists for Q2N if N is odd and a multiple of 3. 
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TABLE I: Qq x Z4 x R assignment of the chiral supermultiplets, where R is the R parity. This is 
an alternative assignment to the one given in [2^. The abeUan Z4 is also an alterative to the abelian 
discrete R symmetry Z12R 26], for which one needs more singlets to construct a desired Higgs sector. 
The anomalies, Qe[SU{2)L]^,Qe[SU{3)c]^, Zi[SU{2)L]'^ and Zi[SU{3)c?, can be cancelled by the Green- 
Schwarz mechanism if, for instance, K2 = K3 is satisfied [ssj. where K2 and K3 are the Kac- Moody levels for 
SU{2)l and SU{3)c, respectively. 



where 



Wc 



Wr 



J2 {yi;'QiU'rH^ + Yf/QiD^^Hf) 

I,i,j,k=l,2;S 



I,i,j,k=l,2,3 

i2 I rTi2 



Xr 



(8) 

(9) 



i=l,2 



Wh = rriTiT't + T^) + myF^ + Xs{Sl + Sl)Y 

+ A3 [-(i/r//f - i^2"^2)ri + {H^H^ + //2"^l')7^2] . 

The Yukawa matrices F's are given by 

/00 0\ /OO K"^'') \ 



Y"3((i3) 



f;^'') 
\o r^"^'^^ y 
/ 



Y"2(d2) 
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(10) 



(11) 
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All the parameters appearing above are real, because we assume a spontaneous CP violation 
to occur. We will shortly come back to this issue. 

The Z4 charges of the fields are so chosen that the most general Qq x Z4 x R invariant 
Higgs superpotential ()10|) has an accidental symmetry: 



where H^''^, T2 and Y do not transform. This symmetry ensures the stability of the VEV 
structure 



(S,) = {S2) = v'e'^''/'\ {T,) = t;^e^(^"/2), (Y) = t;^e^(^"/2), {T,) = 0, (15) 

where f 's are non-vanishing real quantities. 

The total Higgs potential consists of the supersymmetric part {D terms and F terms) 
and the Qq x Z4 x R invariant SSB part, where the B and A terms are given by 

^^SB = Brif^ + f|) + ByY' + XsAsiS', + §1)? 



^3 



+A3A3 



-{H^Hf - H^Hi)Ti + {H^H^ + H^Hf)T2 + h.c. (16) 



(The fields with a tilde are the bosonic components of the corresponding supermultiplets.) 
As in the case of the supersymmetric part, 5's and A's are assumed to be real. We have 
investigated the minimization conditions for the angles ^'s given in ()15|). and found that a 
nontrivial solution of the conditions can exist. How the mass spectrum for this nontrivial 
solution looks like is still an open problem. A complete analysis of this problem, which is 
similar to that of [s^ , will go beyond the scope of the present paper. We will publish the 
analysis elsewhere. To proceed, we here simply assume that there exist a nontrivial CP 
violating set of VEV's. 

C. Fermion mass matrices and diagonalization 

We assume that VEVs take the form (jl5|) . from which we obtain the fermion mass 
matrices. 
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1. Quark sector 



The quark mass matrices are given by 



m 



m 



c f 3 e '"3 



\ -Y^vle-^^^ Y^vle-^^^ ^Yy^e 
( 





^/TY^vie-'^i 







I 



Yifvfje- 
YMe~ 



\ 



(17) 



(18) 



y^Y^vie-'''^^ ) 



In the present case, the unitary matrices that rotate quarks have the following form: 
^l(r)Pl(^ii^0^j^'^, where O's are orthogonal matrices, and 



a. 



uL{R) 



1 
V2 



Rl = 



Rr = 






V2 
\ 



^R 



/ 1 1 \ / 

-1 1 

^ V2 y 

/l 

exp(z2A0") 
^0 exp(zA^") ) 

I exp(i2Ar) \ 

1 
^00 exp(zA^") y 



1 -1 \ 
1 1 
y y 



exp(z^3^) 



and similarly for the down sector. The phase matrices -P^'^ can rotate away the 
the mass matrices, so that we can bring m" into a real form 

/ qjvu \ 

m« = pfRjm''RjiP^ = mt 



\ 





-Qu/yu 
b' 



(19) 
(20) 

(21) 

(22) 

phases of 

(23) 



ylj 



which can be then diagonalized as 



Ofm^Ol 



and similarly for m'^. 

The CKM matrix Vckm is given by 



^ TTlu ^ 

mc 
V mt J 



(24) 



CKM 



OfPl^PtOt 



(25) 



The form of the mass matrix is known as the next-neighbor interaction form | 40l |4 
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and for the set of the parameters 



= ej-O'^ - 61 + r = -1.25, qu = 0.0002143, 6„ = 0.04443, b'^ = 0.09338, 
yu = 0.99732, = 0.005091, ba = 0.02570, b'^ = 0.77606, yd = 0.7940, 

we obtain 



(26) 



m^/nit = 1.10 X IQ-^, mjmt = 4.16 x 10"^ md/m;, = 1.22 x 10"^ m^/mb = 2.13 x 10" 



0.9747 0.2236 0.0040 \ 



\V( 



CKM 



0.2234 0.9738 0.0421 
0.0093 0.0413 0.9991 / 

The experimental values to be compared are 



sin2/3(0i) = 0.738. 



(27) 



42 



iT^exp I 



/ 0.9739 to 0.9751 0.221 to 0.227 0.0029 to 0.0045 \ 

0.221 to 0.227 0.9730 to 0.9744 0.039 to 0.044 
^ 0.0048 to 0.014 0.037 to 0.043 0.9990 to 0.9992 y 



sin2/3(0i) = 0.736 ±0.049. 
The quark masses at Mz are given by 



(28) 



mu/rrid = 0.541 ± 0.086 (0.61) , mjmd = 18.9 ± 1.6 (17.5), 

= 0.73 ± 0.17 (0.72) GeV , = 0.058 ± 0.015 (0.062) GeV, 
mt = 175 ± 6 GeV , m;, = 2.91 ± 0.07 GeV, 



(29) 



where the values in the parentheses are the theoretical values obtained from (jTTj) for rrit = 174 
GeV and mi, = 2.9 GeV. So, we see that the model can well reproduce the experimentally 
measured parameters. Because of the family symmetry, the CKM parameters and the quark 
masses are related. In Fig. 1 we plot the predicted area in the sin 20i — 03 plane. We see 



72 



^^^^^^ 
JSaS^ <^ 



0.6 0.62 0.64 0.66 0.6S 0.7 0.72 0.74 0.76 0.78 0.8 



sin 2c^j(P) 

FIG. 1: Predicted area in the sin2(/)i — 03 plane. The vertical and horizontal lines correspond to the 
experimental values, sin20i(/3) = 0.726 ± 0.037 and (ps = (60° ± 14°) QQ- 
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from Fig. 1 that the model requires 



> 68". 



(30) 



Another quantity to be compared may be |Vtd/l^<(|, whose experimental value has been 
recently obtained from the observation of 6 — > + 7 in the B decays |45|: 



Model : \Vtd/Vts\ = 0.21 - 0.23, 

+0.026 



Exp. : \Vtd/Vt 



tsi - 0.200 (exp.) (theo.). 

' -0.025 -0.029 ^ 



(31) 



We may conclude that 9 independent parameters of the model can well describe 10 physical 
observables. 

The mass matrices given in the present model can be analytically diagonalized and 
the approximate formula for Vckm implies that 



V,,,. 




cb 
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0.81(^)e 

TUb 



nir. 



cd 



— + 0.794W— e 

rrir V m 




■il.25 



(32) 
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which should be compared with the Fritzsch formulas 

Finally, we give the unitary matrices that rotate the quarks for the choice of the param- 
eters given in (pUj) : 



UuLiR) - R1(R)Pl(R)0l{R)^ UdL{R) 

where and ^ are given in (fT^ and (^01) , and 



R 



L{R)^L{R)^L{R)i 



(33) 



' ^R — 



\ 



J 



/ -0.9987 0.0515 -10"^ \ 
0.0513 0.9943 0.0933 
-0.0048 -0.0932 0.9956 

-0.9552 0.2960 -0.0001 \ 
0.1866 0.6025 0.7760 
-0.2298 -0.7412 0.6308 y 

Using these orthogonal matrices and the matrices defined in (fT^ ~ , Eq. gives the 
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(35) 
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unitary matrices in the explicit form: 

/ 0.706 0.0363 1.42 x 10-^ 
-0.706 -0.0363 -1.42 x 10"^ 




uL 



\ 



( 



+ e 



2iAe" 



0.0313 
0.0313 



+ e 



0.0363 -0.705 
0.0363 -0.705 
y -0.00229 e-^^^" 0.0442 e"^^^" 0.999 e''^^^ J 

( 0.695 0.128 0.00352 \ 
0.695 -0.128 -0.00352 


/ 0.128 
0.128 
^ -0.00783 e 

/ 0.0363 0.703 0.0660 \ 



J- 



-0.695 
-0.695 
0.0150 e-'^^" 



0.0115 
0.0115 



1.00 e 



UuR = 



+ e 



UdR — e 3 



0.0363 0.703 0.0660 


/ -0.706 
0.706 



J 



+ e 



y -0.00482 e 

/ 0.132 0.426 0.549 
0.132 0.426 0.549 


/ -0.675 
0.675 



0.0364 
-0.0364 
-0.0932 e-^^^" 



-6.72 X 10-6 
6.72 X 10-6 
0.996 e-'^^" 



v 



2iAe''-+iei 



0.209 
-0.209 



-7.35 X 10-5 
7.35 X 10-5 



\ -0.230 e-^^^' -0.741 e-'^^' 0.631 e-'^^' 



The unitary matrices above will be used when discussing the SUSY flavor problem 
III and proton decay in sect. IV. 



2. Lepton sector 



The mass matrices in the lepton sectors are: 



m 







d „-i6 
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c c 

Yi^ Y'^ 

c c 

\ Y^ Y^ V2yj^tan7"e^(^3- 
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where tan 7" = v^/v'^. We start with the mass matrix of the charged leptons m*^: 

/ \ 







(42) 



y rur J 

One finds [? ] that UeL and UeR can be approximately written as 



U, 



eL 



eR 



R 



R 



I -ee(l + el) -(l/v^)(l elel) 1/v^ \ 

ee(l-e^) {l/V2){l-el + ey^ 1/v^ 
1 - el -V2ee V2eeel J 



V 

/ 1 - ^^/2 



\ 



-4(1- el/2) 







V 



(43) 



(44) 



ele, 1 - el/2 J 



where R interchanges the first and second row, 

/O 1 o\ 
R= 10 
{0 ij 



and = nif^/mr and eg 
becomes 



me/{V2m^). In the hmit rrie 



(45) 



0, the unitary matrix UeL 



1 I/V2 l/V2\ 

-1/V2 1/V2 , 

Vi ; 

which is the origin for a maximal mixing of the atmospheric neutrinos. 



47l | takes place. As 



As for the neutrino masses, we assume that a see-saw mechanism 
we can see from ©, the mass matrix for the right-handed neutrinos is diagonal: rriAr = 
((— l)^mAr, (— l)''mAr, Aatw ^ exp(z6'y/2)), (77 = or 1), where mj^ and \nv^ are real and 
positive by assumption. Therefore, the Majorana mass matrix for the left-handed neutrinos 
becomes 



m (m 



/2(p2)' 



2p2P4 

2(p2)' 
y 2p2P4 2(p4)2 + (p3)2expi2^3 j 



R. 



(46) 



where R is given in ()45|) . and 



2v93 = 2(e^ - r) + r/TT - ^^/2 



P2 = iF;i;]^,/v/^ , P4 = \y/:,vI/^, P3 = -^r:t;3Vv^ 



Y 



(47) 
(48) 



11 



The factor (— l)''e*^^" has no effect, and so we ignore it in the following discussions. 

Noticing that the p's in (I46p are real numbers, we recall that that Mj, can be diagonalized 

as 



m 



\ 



mj^2 



(49) 



\ y 



where Ci2 = cos6'i2, and Si2 = sin 6*12, and 

= R -1 

^ 0126"**^''^'""'"*'^^^/^ 5-^26"*^''^''+*'^^^''^ y 
mi,3 sin = sin ip2 = m^,^ sin ipi , 2^93 = v^i + "^2, 

1 



Am23 sin^ 26'i2 cos^ 0i 
It turns out that only an inverted mass spectrum 



tan (f)y for |r| << 1. 



(50) 

(51) 
(52) 



(53) 



is consistent with the experimental constraint |Am2i| < |Am23| in the present model. Note 
that Eq. ()5H1 is satisfied for 



(54) 



and not for cpi ^ (p2. 

The mixing matrix Vmns can be obtained from UliJJy, where UeL and Uu are given in 
and (jnOl), respectively. The product Ul^Uy can be brought by an appropriate phase 
transformation to a popular form, which in the present model approximately assumes the 
form 



(55) 
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-S12IV2 
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with 
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Sl3 



3.4- ■■ X lO"'', (5, 



CP 



\/2m^ 

sin2Q; = sin(y9i — ^92), sin 2/3 = sin((y9i — ip 



(56) 
(57) 



where ^1,^2 and are defined in There are seven independent parameters to 

describe 12 parameters (3 + 3 = 6 masses, three angles and three phases) of the lepton 



^ Unfortunately, this value of S13 is too small to be measured |4l 
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sector. Therefore, the effective Majorana mass < > in neutrinoless double (3 decay, 
for instance, can be predicted. In Fig. 2 we plot < rriee > as a function of sin^j, for 
sin^^ia = 0.3, Am^^ = 6.9 x IQ-^ eV^ and Am^g = 1.4,2.3,3.0 x IQ-^ eV^ 0. As we see 
from Fig. 2, the prediction is consistent with recent experiments 

We will use the unitary matrices given in (pSjl . ()44p and (jSUjl in sect. Ill and IV, and will 
see that the small parameter S13 = = me/(v^)/m^ ~ 0.0034 appears as a suppression 
factor in FCNCs as well as in some of proton decay modes. 
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^ 0.06 
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0.01 


I 

FIG. 2: The effective Majorana mass < > as a function of siiKp,^ with sin^ ^12 = 0.3 and 
Am^i = 6.9 X 10~^ eV^. The dashed, soUd and dot-dashed hnes stand for Am23 ~ 1-4,2.3 and 
3.0 X 10"'^ eV'^, respectively. The Am^x dependence is very small. 
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51^. 




III. THE SUSY FLAVOR PROBLEM 



If three generations of a family is put into a one-dimensional and two-dimensional irreps of 
any dihedral group, then the soft scalar mass matrix for the sfermions has always a diagonal 
form: 

/ a% \ 



m 



{q/)LL 



mi „- 



\ 










b' 



\ 





q/ 



™2 



m 



uRR 



mi „- 

q,l 







u, d, e) 



(58) 



V 6?,/ 



where m-i denote the average of the squark and slepton masses, respectively, and 
{o.L{R),bL(^ji)) are dimensionless free real parameters of 0(1). Further, since the trilinear 
interactions {A terms) are also Qe invariant, the left-right mass matrices assume the form 



{^Ilr)^. = ^",(m"),, {a = u,d,e), 



(59) 
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where A?'s are free parameters of dimension one, and the fermion masses m's are given in 
(fTTj) . (fT5|) and (^01) . Here we assume that A^'s are in the same order as the gaugino masses. 
They are real, because we impose CP invariance at the Lagrangian leveL 
We work in the super CKM basis and calculate 

= ^l,i)LL UaL and A^^(^^) = f/f^^^^ ALfl(LiJ) UaR (60) 

to parameterize FCNCs and CP violations coming from the SSB sector, where the unitary 
matrices f/'s are given in (j^Hj) . (pSj) and In doing so, one observes that something 

interesting happens; a phase alignment. This is because the only source for CP phases 
comes from VEVs To see the phase alignment, we first observe that the matrices Rl,r 
and the phase rotation matrices P^'r, given in (|T9|l . (j20|) and (j2T|l . commute with the scalar 
soft mass matrices m^/^^ and m^ rf^/j. This implies that A^'^^^ are real, where ^1l,rr 
trivially real as we see from ()43|) and (jUj). As for the left-right soft mass squared we 
find that (-P£)^ {RD'^ih^ij^ RrPr is a real matrix for all a = u,d,e. Consequently, no CP 
violating processes induced by the SSB terms are possible in this model, satisfying the most 
stringent experimental constraint coming from the EDM of the neutron and the electron 



In [101, experimental bounds on the dimensionless quantities 

^LL,RR,LR = ^LL,RR,LR/'"T'le {a = U,d,e), (61) 

are given. The theoretical values of (5's for the present model are calculated below, where 
Aaf = af - 6f , Aa^ = a^- b% = ^ (a = rf, e), (62) 

are introduced. 



Leptonic sector {LL and RR): 



(^12: 


)ll 


= {^21, 


)ll 




4.9 X 10"^ 


Aal, 


(^13: 


)ll 


= (^31; 


)ll 




-1.7 X 10" 


' Aal, 




)ll 


= (^32! 


)ll 


rv-' 


8.4 X 10"* 


Aal, 


(^12: 


^RR 


= ('^21; 


)rr 




8.4 X 10"* 


Aa^, 


(^13: 


^RR 


= ('^31; 


)rr 




5.9 X 10"^ 


Aa^, 


(^23; 


^RR 


= (^32! 


)rr 




-1.4 X 10" 
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Leptonic sector (LR): 

KU.. 2.6X10- (Al-At,) (nifX) 

W3).. ^ 3.1 X 10-' (is, -is) (nifX) 

(523)i:R ^ -1-5 X 10-9 (i^, - Al) ( 
^ -2.5 X 10-« (i^ - At) 



Up quEirk sector {LL and 



(5^2; 


III = 


iS-2l) 


III 




1.0 X 10-^ 


Aai, 




III = 


('^31) 


III 




2.3 X 10-3 


Aai, 


(^23; 


III = 


(^32) 


III 




-4.4 X 10" 


-2 Aai, 




RR = 




^RR 




-4.5 X 10" 


Aa^, 




RR = 


(^31) 


^RR 




4.8 X 10-3 


Aa^, 


(<^23) 


RR — 


(<^32) 


^RR 




9.3 X 10-2 


Aa^. 
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Up quairk sector (LR): 



i^33)LR 
{^12)lR 
iS'i3)LR 
iSlil)LR 
{^23) LR 

{S^,)Ln ~ 10 



10"^ 
10-3 



3.8 (a: -Al- Al + 2Ai) 



500 GeV\ 



-) 



-1.4 [aI -Al- A",) + 7.7 X lO-^i^f ] (^^L^€\^ 



O.SSA" + 6.8 X 10-M." + 3.0 X IO^^A;:, + 1.7 x lO^^M^ 



500 GeV 



■nin 



-{^l^LR ^ 7.4 X 10-5 [K-K- K' + K 



500 GeV\ 



V' 



10-^ 
10-3 
10-^ 



7.9 - - 7.0 X 10 



{aI - Al) 



500 GeV 



-1.7 (i-i^,)- 3.3 X 10-3 (^), 
1.5 [aI - i^) + 1.3 X lO-'i^, + 3.6 X lO-^i^] {^^^^^ 



-3.2 [Al-AIA- 6.4 X 10-3i" - 1.7 x lO-^i" 



500 GeV^ 



Down queirk sector {LL and RRu): 

iSt2)LL = iSi)LL ^ 1.2 X 10-^ Aal, 

iSt3)LL = {Si^)LL ^ 7.8 X 10-3 Aal, 

(Si)LL = (5'32)ll ^ -1.5 X 10-2 Aal, 

{Sf2)RR = {S^i)rr - -1.7 X 10-1 Aai; 

iSn)RR = iSi)RR ^ 1-4 X 10"^ Aa^, 

{Si)RR = {S^32)rr - 4.7 X 10-1 Aaj,. 
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Down quark sector (LR): 



{5i 



llJLR 



{6' 



22)LR 



33 )LR 



LR 



21) LR 



{5: 

iS'^3)LR 
(^31) LR 



{6' 



23 )LR 



32) LR 



10"^ 
10"^ 
10"^ 



6.6 At -At- At + 2A 



500 GeV 



A.l{At-At,) +llAt + 1.5At 



nig J 
■500 GeV 
rrin 



10 



-3 



2.3if + 1.5 X lO'^i^ + 3.5At, + 1.4 x lO'^Ai 

500 GeV\ 



2.1 At-At-At, + A^, 



1.3 (At-At^+SAiAt-A^ 



nig J 



-1.8 {At- At] -2.8 i At, - At 



rrin 



500 GeV \ 
Trig J 
500 GeV 



500 GeV\ 
nig J ' 



-8.4 [At-At,] - 10-^ ( 5.6^4;; - 6.1v4^ 



500 GeV 



nig 



3.5At - 9.3At + 5.3i^, + 0.52if 



500 GeV \ 

nig J 



-2.7 [At-Ai]- 1.8 X 10"^if - 6.7 x lO^^i^^ 



500 GeV^ 



nig 



(68) 





Exp. bound 


Model 


^J\MS^2)iL,RR\ 


4.0 X 10"^ rhq 


(LL)1.2 X 10-^Aa^,(i?i?)1.7 x IQ-^Aa^ 


^\Rei6f,)LLi6f2)RR\ 


2.8 X 10^3 rhq 


4.5 X IQ-''^ ^ AalAaj^ 


\/\MSf2)lR\ 


4.4 X 10~^ rhq 


~ 2 X 10-5 (i^ -At- At, + if)m^^ 


\/\MSf3)lL,RR\ 


9.8 X 10-2 rhq 


(LL)7.8 X lO-^Aa^, {RR)1.4. x lO-^Aa^ 


^J\Rei6f,)LL{6f,)RR\ 


1.8 X 10^2 rhq 


3.4 X lO-^yAa^Aa^ 




3.3 X 10^2 rhq 


~ 2 X io-^{Ai -At + At, - At)m^^ 


^J\MSr2)iL,RR\ 


1.0 X 10~i rhq 


(LL)l.O X W-^Aal,{RR)'i.5 x lO-'^Aa^ 


y^\Re{6^2)LL{Sf2)RR\ 


1.7 X 10-2 rhq 


2.1 X lO-^^^Aa^Aa^ 


\/\MSh)lR\ 


3.1 X 10^2 ^_ 


~ 7 X 10-5 (i;^ -A^- Ai, + i^)mri 


1(^23) ll,rr\ 


8.2 m? 


(LL)1.5 X 10-2Aa^,(i?i2)4.7 x lO-^Aa^ 


\iSl3)LR\ 


1.6 X 10-2 m? 


~ 5 X 10-5 (if - ^fe + At, + O.lif )mri 



TABLE II: Experimental bounds on (^'s and their theoretical values in Qe model, where the pa- 
rameter "rhq denotes mq/500 GeV, and Aa^^/j and A are given in 1)62^ . 



In Table II and III, theoretical values of certain 5's calculated above and their experi- 
mental bounds are summarized. We see from the tables that to satisfy the experimental 
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Exp. bound 


Qq Model 


m2)LL\ 


4.0 X 10"^ m? 


4.9 X 10-3Aa[ 


\{^12) Rr\ 


9 X 10-^ m? 


8.4 X lO-^Aafj 


m2)LR\ 


8.4 X lO"'^ m? 


5.1 X 10-*^(if, - A^^)mJ^ 


I(^!3)ll| 


2 X 10^2 ^2 


1.7 X lO-^Aa^ 


\{^iz)rr\ 


3 X 10^1 m? 


5.9 X 10-2 Aa|j 


m3)LR\ 


1.7 X 10-2 m| 


3.1 X 10-\Al - Al,)mJ^ 


mz)LL\ 


2 X 10^2 m| 


8.4 X 10-^Aa£ 


\{^2'i)RR\ 


3 X 10"^ m? 


1.4 X lO-^Aafj 


\{^2?,)lr\ 


1 X 10^2 ^2 


1.5 X 10-^{Al - Al,)mJ^ 


\{Sh)LL{5h)LL\ 


1 X 10~^ m| 


1.4 X 10-12 (Aai)2 


\{^h)RR{^h)RR\ 


9 X 10"^ m| 


8.4 X 10-8(Aa|j)2 


m:,)LL{5h)RR\ 


2 X 10"^ m? 


5.0 X lO-^Aa^Aa^ 


mz)RR{5h)LL\ 


2 X 10"^ fn) 


2.4 X 10-"Aa[Aa|j 



TABLE III: Experimental bounds on 6^s and the theoretical values in Qq model, where the param- 
eter fh^ denote m^/100 GeV and AaL,R and A are given in (|62|) . 

constraints, the SSB parameters Aa^, Aa£, Al, and of the present Qq model should 
satisfy 

Aa^ < 10-\ Aai < 10~\ Al, - A^ < 10'^ (69) 

while the other SSB parameters are allowed to be of 0(1). So, one can fairly say that Qq 
symmetry can soften the SUSY flavor problem. Note also that the degree of degeneracy of 
the left-handed squark masses Aa^ does not need to be very accurate. We find that because 
of the constraint on Aa^ given in (j69|) . Aa|^ < 10+^ is sufficient to satisfy all the constraints. 
This has an important consequence for proton decay as we will see in the next section. 

IV. PROTON DECAY 

In this section we consider proton decay, which is a process reflecting the flavor structure 
of a model. 

A. Qq invariant baryon and lepton number violating operators 

In supersymmetric models, the baryon and/or lepton number violating operators of lower 
dimensions are l27l: 
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1. dimension-four R parity violating operators, and 



2. dimension-five baryon and lepton number violating operators. 



Both operators are controlled by the flavor structure of a model. As for the dimension-four 
R parity violating operators, Qq flavor symmetry in the model considered in the previous 
sections allows only lepton number violating operators: 



LLE"" 
LQD^ 
H'^NH'' 
NNN 



^iisLiLjE^ + X^iiL^LjEj + X^^^^L^L^E^ + XfuxLjLjE'j^, 

XswH^^NjHf + Asss-f^lA^sifg + XfuKHfNjH^, 
~X'fjjKNjNjNK, 



(70) 



and those in which Lj is replaced by Hf. Note that the baryon number violating opera- 
tor U'^D'^D'^ is forbidden by symmetry. Therefore, the dimension-four operators in the 
present model can not mediate proton decay. 

We next look at dimension-five operators. The baryon number violating dimension- 
: ive operators (which are allowed by R parity and B — L symmetry) can be written as 



27 



28 



5iEi 



1 

M 



E 

j j,fc,i=l~3 



Cf'Q^Q.QuU + CT'U^EpilD 



ij Ml 



(71) 



where the first and the second term is called the LLLL and RRRR operator, respectively. 
In grand unified theories (GUTs), effective dimension-five operators can be generated by 



integrating out colored Higgs multiplets 



27 



and therefore the size of the coefficients of 



the operators strongly depends the Yukawa matrices. For the minimal SUSY SU{5) GUT 
one obtains 



M = Mhc, C 



ijkl 



Cf' = Yij^Y^' at Mgc/t, 



(72) 



where Mhq is the colored Higgs mass of order of the GUT scale, and Yjj, Yd are Yukawa 
coupling matrices appearing in the superpotential 



Y 



^YfflO'lO^H + V2Yi^lO'5^H. 



(73) 



Unfortunately, the minimal SUSY SU{5) GUT should be excluded by the decay mode 
p — > K^u if the gauge coupling unification should be strictly satisfied j29|. 

If we do not assume any GUTs, the baryon number violating operators could be generated 
by some unknown Planck scale physics. In this case, the mass parameter M in fl7ip is given 

by 



M = MpL = 2 X lO^^GeV 



(74) 
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while the coefficients C^"*^ remain undetermined. So, the operators are supphed with a 
suppression factor 1/Mpi which should be compared with I/Mhq — 1/10^^ GeV~^ in the 
GUT case. However, this suppression is not sufficient to keep the proton stable, unless that 
the coefficients C's are smaller than ~ Of 10 "J). An efficient tool to suppress or to forbid 
proton decay is symmetry. Many authors jsS, 30, 5^ have proposed a model in which both 



this is 



the gross structure of the baryon and lepton number violating operators and the structure 
of the Yukawa couplings are fixed by a single flavor symmetry and its breakin g. I 
realized, the flavor symmetry can be tested by proton decay, too. In the models of |23l . 
a certain set of flavon fields is needed to form invariants, and the flavor symmetry is assumed 
to be broken at a superhigh energy scale. In contrast to these models, Qq flavor symmetry 
is broken (spontaneously and at most softly) only at a low energy scale which is comparable 
with the SUSY breaking scale. So, it is natural to assume that flavor symmetry is 
intact at the Planck scale, too; we do not have to introduce flavon fields. Moreover, Qq is 
nonabelian, which can indeed reduce the number of independent coefficients drastically, as 
we will see now. We find that the relevant superpotential containing the baryon number 
violating invariant dimension-five operators generated at the Planck scale can be written 
as 

/,J=1,2 

where the superfields in ()75|) are in the flavor eigenstates. To obtain Eq. lfTSj) . we have not 
assumed any symmetries such as R-parity and B — L except for Qq symmetry, where the 
Qe assignment is given in Table I. 



B. Gross structure of the dimension-five operators and the lowest order approxi- 
mation 

As we can see from flT^I) . allows only three independent coefficients, Cl,C'^^ and 
C'^'. We will see moreover that the ffist term, the LLLL operator, gives the most dominant 
contribution to proton decay, while the RRRR operators can be neglected in the lowest order 
approximation. Consequently, the relative size of all the partial decay rates is fixed in this 
approximation, once the SSB sector is fixed. To begin with, we recall two basic facts: 

1. Since the operators in the ffist two terms in W^^ contain quark fields of the third 
generation in the flavor eigenstate, small mixing parameters appear when fields are 
rewritten in terms of the mass eigenstates, that is, 

^i = V^i^^j\ 1 = 1,2, (76) 

where the subscripts / and m denote the flavor and mass eigenstates, respectively. 
The mixing parameters V37 will be multiplied with the coefficients C's in the decay 
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amplitudes, so that the condition on C's can be relaxed if V31 are small (see (IHU jl ^ fjH^ ). 
Note that within the framework of the MSSM, there is no such suppression. 



2. In most of models, the universality of the SSB parameters at the GUT or Planck scale 
is assumed to suppress the SUSY contributions to FCNCs. This assumption has an 
important consequence that the gluino contributions to proton decay are negligibly 



small if all the squark masses are degenerate 55|. (In Fig. El we show a gaugino 
dressing diagram which contributes to an effective four-fermion operator.) However, 
in our case, we do not assume the universality. We have seen in the previous section 
that the degeneracy of the squark masses of the first two generations is almost exact 
due to Qq symmetry. Moreover, to suppress the SUSY contributions to FCNCs the 
degeneracy of the first two and third generations needs not to be accurate (see Eq. 

and the discussions below). This means that the gluino dressing diagrams may 
not be negligible in our case; we will have to investigate it. 

We now argue that the first term, the LLLL operator, is the most dominant one. It 
is known that for the RRRR operators, the dominant contributions come from the gluino 
dressing diagrams. From the superpotential (f?^ . we first observe that the second term in 
fl75|) . the first RRRR operator, contains two quark fields of the third generation, implying 
that two small mixing parameters will be multiplied when going to the mass eigenstates. 
Further, the third term, the second RRRR term, contains only quark fields of the first two 
generation. That is, the gluino dressing contributions vanish because the squark masses of 
the first two generations are almost degenerate thanks to Qq symmetry. Thus, the LLLL 
operator is the only one which should be considered in the lowest order approximation, as 
we will do it in the following discussions. 

Note that the LLLL operator contains two third generation fields Q3 and L3, but the 
(3, 1) element of the mixing matrix Uei (see fl43|) ) is equal to one, so that it does not act as 
a suppression factor. The dominant diagrams for the LLLL operator are those with gluino 
dressing. The zino and photino dressing diagrams have the same structure as the gluino 
ones, but they are negligibly small because the corresponding gauge couplings are small. So 
we ignore them in our calculations. We have calculated the higgsino dressing contributions 
and found that they can be neglected, too, if tan/? < 10. So we assume this to simplify our 
calculations. We may further approximate that the squark masses are diagonal in the super 
CKM basis. To see this, we recall that the squark mass squared in the super CKM basis 
can be written as 

^a? h « = «'^'e (77) 

LR RR J 

where is the diagonal fermion mass matrix of the fiavor a, and A's are given in (j6U|) . 
As we can see from (jM|l . ()66|) and the nondiagonal elements are sufficiently small. 



21 



FIG. 3: One-loop diagrams contributing to the the effective four-fermi Lagrangian (|78j) . We 
consider only gluino (g) and wino {w) dressings in the lowest order approximation. 

The effective four-fermi Lagrangian C^q can be obtained from the diagrams shown in 
Fig. 121 , where as argued we consider only gluino and wino dressings. We find ^ 



M^d,s M,N^d,s 
l — e,fj, l — e,fj,,T 



CLL^uduef^^' = Cl^^uduef^^' + CfL^uduef^^', 
CLLiudduY^'''' = Cl^iuddvf'''''' + CtL{uddvY'^''\ (79) 

where and stand for the contributions coming from the gluino and wino dressing 
diagrams, respectively. These are explicitly calculated to be 

CUuddvY'^^^ = 47r^ Yl UHU^lUI^UH {F^{al,bl) - F'^{a%al)) , 

/=1,2 

CUuddvY'^''^ = Ana, ^ [UiM^Ul^U!l{F^{alal) + F^iblX)) , 

7=1,2 

- U'^U^m^U^l {F^ial hi) + F^{a% fel))] , (80) 



CUuduer^^ = -An^ ^ ^IMlU'^IuH {F~^{a%hl) - F^(ai,ai)) , 

7=1,2 

CUuduef'^'^ = -Ana, ^ MMl U^lu!^ {F^ (al al) + F^ (bl b^)) , 

7=1,2 

- UiMtUilUH {F^ial bl) + F^ial, b',))] , (81) 



^ We use the notation of l57i IS; 
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where the loop functions are defined as j57|, |5 



with 



1 1 



Xgi Inxgi _ Xg3 lnxg3 



Xwi 1 



m 



m 



q q 

'g,w 



9 

mi 



m 



2 — 

g,w 



2 

m- 



m 



We also use Ag = Aa^ 



(which is first introduced in (jU^ ) and 



(82) 



(83) 



Then we calculate the decay amplitudes as a function of Xg^ for given values of rgAg and 
7/^3, while for simplicity we assume the GUT relation between the wino and gluino mass ® 



rriu, = 0.27mg. 



(84) 



The unitary matrices U(u,d,e)L in C'll are explicitly given in (jHUj) . (jHTj) and PHj). where the 
individual phases ^"''^ and ^3'^^ in the quark sector are not fixed fsee()26|)). However, it is found 
from (jHEI) and (jSH) that the phase dependence of the combinations appearing in Cll, that 
is, Ull^Uli^''^^ and U^^^^^U^^^^j^, is small, and moreover the absolute size of the suppression 
factor U^^j^^j^ is independent of the phases. Therefore, in the following calculations we choose 



= 91 = 0, ei 



-1.25 



(85) 



without significantly changing the results. 



C. Gluino versus wino contributions 

Before we calculate the decay amplitudes, we investigate the relative size of the gluino 
and wino contributions. As pointed out, to suppress the SUSY contributions to FCNCs in 

model the degeneracy of the first two and third generations needs not to be exact. In 
fact, as we can see from Table III and Eq. (jH^jl . Ag < 10 is sufficient. (Ag = Aa| defined 



c masses.) If this is a case, the 
5^, so that they may dominate 



in (jU^ expresses the degree of the degeneracy of the squar 
cancellation of the gluino contributions is no longer perfect 
over the wino ones. In Figs. lUini we present the results on the ratio & jC"^ as a function 
of Xgz for r^Ag = (10,1,0.1,0.01) and y^z = 10 (right), 1 (left). Figs. IHshow the results 



If we change the GUT relation the results below can change only slightly. 
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on C{uddvys which control the size of the anti-neutrino modes, while Figs. El are those on 
C{udue)^s for the decay amplitudes into a charged lepton. We can see furthermore that 
(except for C{udueY'^^^ shown in the last two figures of Figs. ^ the gluino contributions 
can become of the same order as the wino ones for Xg^ ~ 0(1) if r^Ag ^ 1. We can also see 
that the smaller the value of r^Ag is, the smaller are the gluino contributions, as it should 
be. 

If the degeneracy of the squark masses is assumed, the gluino contributions vanish, and 
the wino contributions are then the dominant ones. This is assumed in most of the previous 
investigations on proton decay. Let us call the wino contribution under the assumption of 
the degeneracy (that is, Ag = 0) Cll{- ■ ■ )^ssmi ^-^^ compare it with Cll{- ■ ■ )gg for which 
Ag may be different from zero. Since it will turn out that the decay mode p K^v has 
the largest rate, we restrict ourselves to Ciiiudduy^'^'- (which gives the largest contribution 
to the decay mode). Fig . Elshows the ratio C^Liiuddv)]^'^^ /Cf^j^iuddvY^^^ggj^j as a function of 
Xg^ for r-gAg = 10 and = 1, 10. We see from the figure that the gluino contribution in Qq 
model is smaller than the wino ones of the MSSM, although the SU{?>)c gauge coupling is 
larger than that the SU{2)l gauge coupling. This is due to the fact that the LLLL operator 
in ()75p contains one which requires a small mixing parameter to appear when expressed 
in terms of the mass eigenstates of the first two generations. 



D. Partial decay widths and their relative size 



Here we calculate the partial decay widths of the proton. The partial decay widths for 
the decay mode p Mi, T{p ^ Mi), can be written as [s^ 



T{p^Mi) = { —^if^ ] '"Z "f'f} \Amp {p ^ Mi)f , (86) 



with 



Amp{p^K+i)) = K^CLdudduf^^^ + K^CLL^udduf^^^ (87) 

Amp{p^'R+v) = V2n^CLL{uddvY^^\ (88) 

Amp{p-^ K^ei) = KiCLLiuduef^^^ (89) 

Amp(p^7r°e+) = k^C LL{uduef^^\ (90) 



where k's are defined 



by 



Ki = l + —{F -D) = 0.70 , K2 = l + — ^(3F + D) = 1.6, 

•niB 3mB 

n, = l + ^(F + D) = 2.0 , K, = 2 + = 2.7, (91) 

rriB TUB 

K, = i=(l + F + D) = 1.Q , K,= '^D = 0.4, 
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FIG. 4: The ratio C'Huddu) / {uddv) is plotted as a function of Xgs for r^Ag = (10, 1, 0.1, 0.01), 
which expresses the degree of the degeneracy of the squark masses, where x^s is defined in (|82|) . 
C^'s stand for the gluino contributions and C""s for the wino ones, respectively. The figures in the 
left column correspond to y^s = 1, and those in the right column to = 10. The corresponding 
amplitudes are responsible for the anti-neutrino modes. 



and 



D 

mp 



0.81, F = 0.44, Pp = 0.003 GeV^ U = 131 MeV, 

938 MeV, = 1150 MeV, ttik = 495 MeV, = 140 MeV. 



(92) 



25 




rgAQ=10 
rgAQ=l 
rgAQ=0.1 
rgAQ=0.01 
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FIG. 5: The same as Fig. |l]for the charged lepton modes. 



0.5 

a 




2 4 6 8 10 



FIG. 6: Ratio of Cf^i^{uddu)^^'' to C]^i^{udduYj^ggj^ as a function of Xgs, where r^Ag = 10 is 
assumed for Cf^i^{uddv)^^'- , where Xgs is defined in (|82j) . C LiiuddvY^"^^ is the largest coefficient 
for p ^ i^+P, and the degeneracy of the squark masses (Aq = 0) is assumed for CYiiudd^Y^ggj^. 



Here, F stand for the coupling constants for the interaction between baryons and mesons, 
I3p for the hadronic matrix element, /vr for the pion decay constant, nip for the proton mass, 
ruB for the averaged baryon mass, and mx and m^r for the Kaon and pion mass, respectively. 
Further, A is the renormalization group(RG) enhancement factor for the coefficient Cl 
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28|,|53|, which in our case becomes 

A = 10.5 (for the MSSM), 12 (for Qg model). (93) 

(Since Qq model contains four more Higgs doublets than in the case of the MSSM, the 
enhancement factor is slightly larger.) 

As the first task, we would like to compare the decay rates in the MSSM and Qq model. By 
the MSSM decay rates we mean the decay rates which are obtained from the LLLL operators 
under the assumption of the degenerate squark masses (no gluino dressing). We also assume 
that all the coefficients of the operators are equal to the single constant Cf^^^^^'^ ^ and that 
all the fields appearing in the operators are in the mass eigenstates (no mixing matrix). In 
Figs. [7J we plot for each decay mode, the experimental bound (thick solid line), the partial 
lifetime calculated in Qg model (solid line) and that in the MSSM (dotted line), where we 
assume 

Ci = Cf^5^ = l, m,- = lTeV. (94) 

The lifetime in model is calculated from the sum of the gluino and wino contributions for 
TgAg = 0.01,0.1, 1, 10, which corresponds to four solid lines. We first see from the figures 
that if TgAq varies from 0.01 to 10, a change of the life time about an order of magnitude 
can appear. We also see that the decay rates in Qg model is much more suppressed than 
those in the MSSM. Quantitatively, we find that the experimental bounds can be satisfied 
if the coefficient for the LLLL operator satisfies 

Cl < 10-(^~5), (95) 

while Cf^issM ^ iQ-id^r) should be satisfied for the case of the MSSM. That is, Qg fiavor 
symmetry can suppress proton decay by four orders of magnitude (which can be also seen 
from the figures). 

Next we would like to compare our results with those obtained in the minimal SUSY 
SU{5) GUT. In our lowest order approximation, only Cl is an independent coefficient, 
implying that the ratio of partial decay widths is independent of Cl- That is, in Qq model, 
the relative size of the partial decay rates is fixed (once the SSB sector is fixed). First we 



recall the case of the minimal SUSY SU{5) GUT 30, 53|. The superpotential for the baryon 



number violating effective dimension-five operators in this case is given by 

^^""^'^ = ^y-y<^^ ^^*^2i (QiQi) {Q2L1) + vcVdi (v*)i, {Q2Q2) {QiLi)] (96) 

in the approximation that the third generation is dropped, where Uu^c-, Udi are the diagonal 
Yukawa couplings of the corresponding quarks, / = 1,2 is generation index (i.e., dl = d 
and d2 = s), and V is the CKM matrix. In writing a nontrivial assumption is made; 
the up quark Yukawa matrix is diagonal over the whole range of energies. Therefore, the 
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FIG. 7: Partial lifetime of the proton for each decay mode as a function of Xg^ with Cl = Cff^^'^ = 
1, nig = 1 TeV. The experimental bound (thick solid line), the partial lifetime calculated in Qq 
model (solid line) and that in the MSSM (dotted line) are plotted. The lifetime in Qq model is 
calculated from the sum of the gluino and wino contributions for VgAg = 0.01,0.1,1,10, which 
corresponds to four solid lines. 

superpotential is not a unique prediction of the the model. Under the assumption of the 
degeneracy of the squark masses, only the wino dressing diagrams contribute to the decay, 
and one finds [sO, 
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T{p-^ K+v) = 7 ( 

r — > TT^Z/) = 7 ( 



2 \2 I 



2\2 



2ycys sin^ 6c tan 9cV2k5 

2 

2 



r -> K^el) = 7 [nip - m\Y lUnUdi (^*)2« cos6'c/ti 
T [p^ -K^ef) = 7 (mp - mlY ly^ydi {V*)2i sin 60^5 

where we have dropped the corresponding loop functions (j82j) . and the common factor 
given by 
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FIG. 8: Ratio of partial decay rates for each decay mode with r^Ag = (10, 1, 0.1, 0.01). Figures 
in the left column are for 1/^3 = 1, and the right column for y^^ = 10. 



From these partial decay widths, we obtain the relative decay widths in the minimal SU{5) 
GUT 00: 
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FIG. 9: Ratio of partial decay rates for each decay mode with r^Ag = (10, 1, 0.1, 0.01). Figures 
in the left column are for 7/^^,3 = 1, and the right column for y^s = 10. 



The corresponding results for Qq model are found to be 
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where as in the case of the minimal SUSY GUT we have suppressed the loop functions. The 
upper (lower) numbers on the right hand side correspond to the wino (gluino) contributions. 
Note that in contrast to the case of the minimal SUSY GUT the mixing parameters explicitly 
appear, reflecting the flavor structure of the present Qq model. The most remarkable dif- 
ference is the charged lepton modes. As we see from ()108|) and ()109|) . they are proportional 
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to |f/e|P, where the unitary matrix Uei (which rotates the left-handed charged leptons) is 
exphcitly given in (j43p . We find 

B{p^ i^^V+) _ B{p^ 7rV+) _ \U^l\^ 
B{p^ K^e+) ~ B{p^ 7rOe+) ~ |f/3i|2 

\ ^ 

— ) ~2.37xl0"^ (110) 
mj 

The same ratio in the case of the minimal SUSY GUT becomes B {p ^ K^fi^) ~ 10'^ x 
B {p ^ K'^e'^) if we use the superpotential (jHEI)- In Figs. |H1 and El we plot four different 
combinations of ratios as a function of Xg^ for various values of TgAq in Qq model. The 
figures on the left (right) side correspond to 2/^3 = 1(10). From the figures, we find the 
hierarchical structure of the decay modes: 

B{p^ ir°e+) < B{p^ 7r°e+) < B {p ^ vr+z/) < B [p ^ K+u) . (Ill) 

In the case of the minimal SUSY GUT we obtain instead 

B{p-^ 7r°e+) <B{p^ A^°e+) « B {p ^ tt+z/) < B [p ^ K+u) . (112) 

Note, however, that although B {p 7r°e+) < B {p —>■ K'^i?), they are basically in the same 
oder in the Qq model. That is, the Qq model predicts that once the decay mode K'^u is 
experimentally observed, then it is likely to observe the decay mode 7r^e~^, too, in sharp 
contrast to the case of the minimal SUSY SU (5) GUT. 




V. CONCLUSION 

Flavor symmetry can play important rolls in supersymmetric models js^. We have 
investigated the SUSY contributions to FCNCs and to proton decay in a supersymmetric 
extension of the SM based on a binary dihedral family group Qq. We have seen that the 
discrete low energy fiavor symmetry Qq can be an alternative to the universality assumption 
of the soft supersymmetry breaking parameters. That is, the existence of a hidden sector, in 
which supersymmetry is assumed to be broken in a fiavor blind manner, is not indispensable 
in this model. Therefore, a variety of supersymmetry breaking mechanisms may become 
phenomenologically viable. 

It has turned out from our analysis on FCNCs that the degeneracy requirement of the soft 
scalar masses of the left-handed squarks in this model can be significantly relaxed. This has 
a considerable effect on the gluino-mediated one-loop amplitudes on proton decay. In most 
of the previous calculations, the degeneracy was assumed so that the gluino contributions 
are cancelled with each other, which implies that the wino contributions are the most dom- 
inant contributions. We have fond that the non-degenerate squarks can change the decay 
rate in the charged lepton modes by an order of magnidude. We have also found that in 
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the lowest order approximation there is only one independent coefficient for dimension-five 
operators that lead to proton decay. Consequently, the relative partial decay rates in this 
approximation are fixed, refiecting the fiavor structure dictated by Qq. 

Our main finding is that Qq fiavor symmetry acts in such a way that the 
smallness of (HiNs)e3, the suppression of /i — > e + 7, and the smallness of 
i? (p — i^°(7r°) + /i"*") /i? ^ i^°(7r°) + e"*") have the same origin; the electron is much 
lighter than the muon and the tau. They in fact vanish in the mg limit. 
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